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RESUMEN 

El resumen sera traducido al espanol por los editores. 

We present a tool, VELNEB_3D, which can be applied to the results of 3D 
photoionization codes to generate emission line profiles, position- velocity maps and 
3D maps in any emission line by assuming an arbitrary velocity field. 

We g ive a few examples, based on our pseudo-3D photoionization code 
NEBIL3D l|Morisset. Stasinska fc Penal feoOS) which show the potentiality and use- 
fulness of our tool. One example shows how complex line profiles can be obtained 
even with a simple expansion law if the nebula is bipolar and the slit slightly off- 
center. Another example shows different ways to produce line profiles that could be 
attributed to a turbulent velocity field while there is no turbulence in the model. 
A third example shows how, in certain circumstances, it is possible to discriminate 
between two very different geometrical structures - here a face-on blister and its 
"spherical impostor" - when using appropriate high resolution spectra. Finally, we 
show how our tool is able to generate 3D maps, similar to the ones that can be 
obtained by observing extended nebulae with integral field units. 

ABSTRACT 

We present a tool, VELNEB_3D, which can be applied to the results of 3D 
photoionization codes to generate emission line profiles, position- velocity maps and 
3D maps in any emission line by assuming an arbitrary velocity field. 

We g ive a few examples , based on our pseudo-3D photoionization code 
NEBIL3D llMorisset et aLl fcoO.^ which show the potentiality and usefulness of our 
tool. One example shows how complex line profiles can be obtained even with a 
simple expansion law if the nebula is bipolar and the slit slightly off-center. Another 
example shows different ways to produce line profiles that could be attributed to a 
turbulent velocity field while there is no turbulence in the model. A third example 
shows how, in certain circumstances, it is possible to discriminate between two very 
different geometrical structures - here a face-on blister and its "spherical impostor" 
- when using appropriate high resolution spectra. Finally, we show how our tool is 
able to generate 3D maps, similar to the ones that can be obtained by observing 
extended nebulae with integral field units. 
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1. INTRODUCTION 

The study of emission line profiles in nebulae 
(planetary nebulae and HII regions) is important 
for many reasons. The shapes of emission lines re- 
sult from the motions along the line of sight of the 
parcels of gas covered by the observational aperture 
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and in which these lines are emitted. With the 
help of some hypothesis, the analysis of line pro- 
files allows one to describe the internal motions in 
nebulae. It also opens access to the third dimen- 
sion in the study of nebulae, since an analysis of 
surface bright ness profiles alone may be mislead- 
ing (see e.g. iMorisset et alJ EoO^I) . A proper de- 
scription of the geometry and kinematics of neb- 
ulae is a prerequisite to understand their dynam- 



1 



2 



MORISSET & STASINSKA 



ical evolution since the epoch of their formation. 
Unfortunately, the interpretation of line profiles is 
extremely difficult. The first spectra that showed 
line splitting in planetary n ebulae were published by 
ICampbell fc Moore] l|l918j) and were interpreted as 
due to rotation. The interpretation in terms of neb- 
ular expansion, which physically makes more sense, 
appeared only in 1931 l|Zanstra|) . An extent discus- 
sion of internal motions i n plan etary nebulae appears 
with the work of lWilsonllll950h . Among other things, 
he noted that the proportion of nebulae with split- 
ted lines in which the blue component is stronger 
than the red component is equal to that of nebulae 
in which the reverse is occurring, which implies that 
the nebulae are transparent in the observed spec- 
tral range. He concluded that the relative strengths 
of the red and blue components depend upon non 
uniform distribution of the nebular material and 
upon the direction of the observer. He also noted 
that the separation between the two components is 
smaller for [Ne V] than for [O III] and [Ne III] and 
concluded that the expansion velocity of the nebular 
material increases with distance to the star. This 
however was shown later n ot to be a general prop- 
erty of planetary nebulae llSab badin fc Hamzaogh; 
Il982j) . bsterbrock. Miller fc Weedmanl l|l966j) ob- 
served emission line profiles in several bright plane- 
tary nebulae, and noted that the lines are broader 
than the expected thermal Doppler widths, indi- 
cati ng significant mass mo tions within the nebu- 
lae. IWilsonl l|l958h and then I Weedmanl l|l968h were, 
to our knowledge, the first to perform a spatio- 
kinematic analysis of the data under the assump- 
tion of ellipsoidal geometry. For the five nebulae he 
studied, Weedman found that the expansion velocity 
increases with distance from the star, and that the 
gradient extrapolates to zero velocity at a point be- 
tween the nebular shell and the center. This implies 
that the expansion velocity has increased since the 
time of ejection. 

Since then, spatio-kinematic studies of plan- 
etary nebulae have gone into several directions. 
One is a purely empirical approach which, from 
observations of splitting, width and intensity 
ratio maps of emission lines infers the nebu- 
lar struc ture, using a few simplifying assump- 
tions l|Sabbadin. Bianchini. Ortolani fc Strafellal 
119851) : this technique has been called to- 
mography. Other studies have focused 
on the kinematics of microstr uctures (e.g. 
Balick. Rugers. Terzian fc Chengalurl Il993l 
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l)2003ft : iGesicki fc Ziilstral l|2003l and references 
therein) have determined the internal veloc- 
ity fields of about 70 planetary nebulae, by 
comparing observed profiles with profiles com- 
puted for ID photoionizat i on m odels, while 
iMorisset. Gruenwald fc Viegasl (|2000^) have devel- 
oped a tool to derive morphological and kinematical 
information on planetary nebulae, based on a 3D 
photoionization code. 

But 3D photoionization codes are heavy in terms 
of computing time, while interpretations of kine- 
matic observational data for real planetary nebu- 
lae with a ID tool can be misleading. Here, we 
present a tool, VELNEB_3D, and its association to 
a quick pseudo -3D photoionization code NEBU_3D 
l|Morisset et al.ll2005|) which allows one to easily ob- 
tain realistic emission line profiles and radial velocity 
maps for any geometry. 

The main goal of the present paper is to empha- 
size a few points: 

• Emission line profiles are a combination of a 
morphology of the emitting gas and a velocity 
field. The complexity of observed line profiles 
can be attributed to a complex velocity law (as 
has been done by Gesicki, Zijlstra et al. for 
years, using simple non-realistic spherical ge- 
ometry) or to a complexe morphology (as com- 
monly observed). 

• Velocity fields can be partially attributed to tur- 
bulence (as claimed by Gesicki & Zijlstra, 2003, 
for a number of planetary nebulae) simply be- 
cause the analysis unduely relies on spherical 
models. If the models account for departure 
from sphericity (as ours do) , there is noo need of 
turbulence to account for the observed profiles 

• Line profiles can allow to distinguish between 
very different morphologies (a blister and its 
spherical impostor) 

The structure of the paper is the following. Sec- 
tion presents the VELNEB_3D tool and compares 
it to previous similar tools. The three next sections 
present a few applications. The main goal of the 
present paper is to present the capacities of this tool 
and no attempt is made to reproduce observations 
of any real object. Section [3] illustrates how complex 
emission line profiles can be obtained using simple 
expansion laws. Section ^demonstrates how aspher- 
ical geometries can, in certain circumstances, lead 
to a broadening of emission lines that can be mis- 
taken for turbulence. Section [5] shows an example 
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when two very different geometries can be distin- 
guished using emission line profiles. Section 6 shows 
an example of a 3D map in several lines generated by 
our code. Section 7 presents a summary and some 
prospects. 

2. VELNEB.3D 
2.1. Description 

VELNEB_3D (which has been developed in IDL) 
uses the results of a 3D photoionization code. In 
this paper, we use the pseudo-3D photoionization 
code, NEBU.3D. Both NEBLL3D and the visuali- 
sation tool VISNEB_3D h a ve be en extensively de- 
scribed by iMorisset et al.l l|2005l) . In this section 
we briefly recall how they operate. NEBU_3D is 
a set of tools that allows the user to construct a 
3D photoionized model nebula from a set of runs 
of a ID photoionizatio n code (in this paper, the ID 
code used is NEBU, se e IMorisset fe PeauignotJll996l 
IPeauienot et al.ll200l[) . The interpolation is done on 
a coordinate cube for all the relevant physical param- 
eters (electron temperature, density and ionic frac- 
tions) and for the emissivities of emission lines the 
user could need. After the construction is done, VIS- 
NEB_3D is used to rotate the nebula to any orien- 
tation, and to compute emission line intensity maps 
integrated along any line of sight. 

One can then attribute a velocity vector to each 
cell of the coordinate cube containing the nebula. 
For example, one can adopt a simple radial expan- 
sion velocity law, where the magnitude of the veloc- 
ity is proportional to the distance to the star, like 
in a Hubble flow. In each cell of the cube, the code 
computes elementary emission line profiles, using the 
velocity vector defined previously and taking into 
account the thermal broadening obtained from the 
local electron temperature and the mass of the emit- 
ting ion. A turbulence velocity can be quadratically 
added to the thermal velocity. An integration of the 
individual emission line profiles along a line of sight 
is performed for each pixel of the surface brightness 
map. The procedures used here are very similar to 
the on es described in more detail in IMorisset et al] 

fcoooh . 

Aperture effects can be simulated by applying a 
mask corresponding to a slit size, orientation and 
shift relative to the center of the image. The effect 
of the seeing is taken into account by smoothing the 
aperture mask by a given angular size square kernel. 
The resulting profile can also be finally convolved by 
a theoretical instrumental profile. 

All this constitutes the VELNEB_3D tool, which 
is here applied to results from a photoionization 



calculation with NEBU_3D, but can be applied 
to results from any photoionization code. As a 
matter of fact, VELNEB_3D in itself is likely very 
simila r to the spatio-kinematic code mentioned 
by lHarman. Brvce. Lopez. Meaburn fc Hollowavl 
<|2004l) . The fact that VELNEB.3D is coupled to a 
photoionization code ensures that the line profiles 
are computed in a consistent way, taking into 
account the distribution of ions and temperature 
throughout the nebula rather than arbitrarily 
assuming a source function for each line. 

2.2. Comparison with similar tools 

We now compare the relative merits and draw- 
backs of several recent tools that have been designed 
to study nebular expansion. 

iFrank. Balick. Icke fc Mellemal {l993) compute 
2D fully hydro simulations including photoioniza- 
tion. Synthetic position- velocity maps are computed 
and compared by eye with real observations, show- 
ing that the 2-wind model accounts qualitatively 
for the diversity of PN shapes and velocity fields. 
This is probably the most elaborate approach to 
study the morpho-dynamics of planetary nebulae. 
ISchonberner. Ja cob. Steffen. Perinotto. Corradi fe Ackerl 
1120051) and ISchonberner. Jacob fc Steffenl ll20f)! ~ 
develop a ID hydromodelling code, including 
photoionization. They compare surface brightness 
distributions and line profiles to observations of 
a few planetary nebulae. One of the big limita- 
tion of these works is th e sphe rical morphology. 
iRiikhorst. Mellema fc Ickel l|2005ll performe three- 
dimensional AMR simulations of point-symmetric 
nebulae and compute PV-diagrams, but only for 
the hydrogen rec ombination lines. In addit ion, both 
this code and the ISchonberner et al code are 

yery d emanding in computing time. Morisset et alJ 
(2000) use a static 3D photoionization code. The 
diffuse radiation treated in the on-the-spot (OTS 
approximation) - similarly to the former two groups. 
They apply their tool to a simple case showing 
how velocity field and geometry go together. The 
main draw-back of this code is the absence of a 
full treatment of the diffuse field, despite the huge 
time of execution needed to comp u te a model. 
lHarman. Brvce. Redman fc Hollowavl l)2003(l used 
a morphological modelling code to help determine 
the geometry (set to axisymmetric), structure 
and kinematics (a linear expansion law is used) 
of the ellipsoidal PN S a 2-21. The code SHAPE 
l)Steffen fe Loped HoOff), which works in 3D, is 
very versatile concerning the explo ration of line 
profiles, but, as for the code of lHarman et alJ 
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(2003), the distribution of the emissivities is ar- 
bitrary and do not come from a self-consistent, 
photoionization approach. It must use the out- 
put of a 3D photoioinization code to improve 
the emissivities distribution. The main draw- 
back of the actual version of SHAPE is the use 
of a commercial 3-D modeling software devel- 
opped to run under the Microsoft operating system. 
ISabbadin. Turatto. Raeazzoni. Cappellaro k Benettil 
l|2006l and references therein) try to solve the 
long-standing problem of "deprojecting the bi- 
dimensional apparent morphology of a three- 
dimensional mass of gas". They use the zero- 
velocity-pixel-column information to constrain 
ID photoionization models and the central-star- 
pixel-line to constrain the expansion velocity 
of the emitting gas. Those technics lead to a 
possible reconstruction of observed nebulae in 3D. 
ISabbadin. Turatto. Cappellaro. Benetti k Ragazzonil 
(2004, and references therein) use tomography to 
reconstruct the morphologies of PNs, but are 
limited by the assumptions of linear expansion laws. 
lOesicki. Ziilstra. Acker. Corny . Gozdziewski k Walshl 
(2006, and references therein for the same group) 
use a ID static photoionization code under the 
OTS assumption to compute emission line profiles 
and compare them to observations of about 100 
planetary nebulae. They allow for complex velocity 
fields, but are severely limited by the spherical 
geometry. The only fully 3D photoionization 
codes with accurate transfer are M OCASSIN 
llErcolano. Barlow. Storey k Liul l2003l) and the 
code of IWood. Mathis k Ercolanol l )2004|h There is 
no computation of line profiles yet, but VISNEB 
can be applied to them 

Ideally, of course, one would like fully hydrody- 
namical, 3D models including photoionization. Even 
when this exists, static pseudo-3D models such as 
ours will remain useful, because of their relative 
simplicity, to explore the parameter space. Of the 
codes explicitly designed for comparison with data, 
the cunjunction of NEBLL3D with VISNEB so far 
appears as the best compromise between relevance, 
versatility and rapidity. 

3. COMPLEX PROFILES WITH A SIMPLE 
EXPANSION LAW 
3.1. A bipolar nebula toy model 
The photoionization model presented in this sec- 
tion consists in applying NEBU_3D to a set of 40 
runs of NEBU, each one corresponding to a differ- 
ent polar angular distance 9 and characterized by a 
different gas density distribution. We consider no az- 
imuthal variation for the density distribution. Since 
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Fig. 1. Electron density distribution, for the bipolar 
model described in Sect 3.1, in which the hydrogen den- 
sity is kept constant along any radial direction. 

the model is axisymmetric, we need to model only 
l/8th of the nebula, obtaining the rest of it by ro- 
tation. Thus, we vary 9 from to 7r/2, from the 
equator to the pole. 

The inner radius, i?inner(0), is defined in such a 
way that the inner cavity has a bipolar shape, de- 
lineated by two equal intersecting spheres of radius 
R2 . The center of each sphere is displaced on the po- 
lar axes by a distance Ri from the ionizing central 
source. The equation for the inner radius is then: 



R\ 



R\ cos 2 (6>) + sin(0) 



We assume that the hydrogen density at the in- 
ner surface is related to the distance of this sur- 
face to the ionizing star by: ni nner (6') = n-cquat x 
(■Rcquat/-Riimcr(#)) 2 - For each polar angle 9, the 
variation of the hydrogen density along the ra- 
dius is given by: n(R,9) R>RiancrW = n inncl (9) x 
(i?inner(0)/-R) 7 (n. b. 7 = corresponds to a den- 
sity that is constant along each direction, 7 = 2 
corresponds to a free expansion in each direction). 
The free parameters for these geometries are then 
: i?i/i?2,-Requat,^cquat and 7. For the model pre- 
sented in this section, we fix R1/R2 so that the 
equatorial inner size of the nebula is half the po- 
lar inner size. This assumption implies 3i?2 = 5R±. 
We also take 7=0, and fix i? e quat and n eqU at so that 
the size of the shell is of the order of the size of 
the inner cavity, namely i? eqU at = 4 x 10 16 cm and 
«oquat = 3 x 10 4 cm~ 3 . 

The stellar energy distribution is that of a black- 
body with T e g =70 kK. We adopt a stellar luminosity 
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of L* = lxl0 37 erg s _1 . The elemental abundances 
are taken solar. Dust is not included in the model. 

The interpolation is done in a cube of 100 3 pixels, 
each of 5xl0 15 cm in size. At an assumed distance 
of 2.7 kpc from the observer, this corresponds to an 
angle of 0.075"for each pixel. 

We show in Fig. ^ a cu t in a plane containing 
the polar axis, with the electron density represented 
by different levels of grey. Black zones correspond to 
null electron density (either because of lack of matter 
like in the inner region, or due to the fact that no 
ionizing photons reach this zone, like in the outer 
parts). From the simple assumptions used to define 
the morphology of the nebula, it follows that the 
shape of the recombination front is bipolar, with a 
strong enhancement of the density at the equatorial 
ring. 

3.2. Computed monochromatic images and line 
profiles 

The images presented in Fig.|5]are computed for a 
nebula with a polar axis making an angle of 45degr 
with the plane of the sky. A seeing of 0.5" is as- 
sumed. The figure shows the surface brightness of 
the nebula in four emission lines: H/3, [O III] 5007A, 
[N II] 6583A and He II 4686A. The bipolar shape is 
seen in all the lines, but the images differ due to the 
ionization structure of the nebula. Obviously, it is 
the [Nil] 6583A image which best suggests the pres- 
ence of the inner cavity. There is a small "noise" 
(the classical "Moire" effect) due to the finite size of 
the cells. This effect is also present in the surface 
brightness profiles shown in the next sections. This 
effect is more important for emission lines that are 
occupying a thin shell than for HI line for example. 
It decreases when the number of cells of the main 
cube increases. 

The emission line profiles are computed for vari- 
ous apertures, indicated in Fig. [2 Solid lines corre- 
spond to centered apertures, dashed lines correspond 
to off-center apertures. The off-center apertures were 
shifted by 0.3" and 0.5" in the x- and y-directions 
respectively. The proportion of the H/3 flux collected 
through the l"x3" and the 2"x2" centered apertures 
are 18 and 19% respectively. 

To compute the line profiles, we adopted a very 
simple expansion velocity law, given by the following 
expression: V(R) = Vq.R/ R max , where R is the po- 
sition vector originating at the centre of the nebula, 
Rmax is maximum distance of the ionization front 
to the centre (i.e. the radius in the polar direction), 
and Vo is set to 40 km s _1 . 

The emission line profiles are computed on a 60 
pixels grid (from -50 to 50 km s -1 ). 



The results are shown in Fig. [3] for three lines: 
H/3, [N II] 6583A and [O III] 5007A (respectively left, 
middle and right panels of the three bottom rows). 
The shape and size of the slit is indicated at the top 
of each panel. Profiles obtained through off-center 
slits are represented by dashed lines. The lowest 
row of panels is for an aperture that covers the en- 
tire nebula. This case corresponds to what would be 
observed for an extragalactic planetary nebula, for 
example. The figure also shows some plots that are 
useful to understand the line profiles. The top panels 
give the radial distribution of the hydrogen density, 
the electron temperature and the expansion velocity. 
Solid curves correspond to the polar direction, dot- 
ted curves correspond to a direction perpendicular 
to it. The second row of panels shows the surface 
brightness distribution in the H/3, [Nil] 6583A and 
[O III] 5007A lines along the polar axis (solid curves) 
and in the perpendicular direction (dotted curves). 

As expected, the lines showing most structure are 
the [Nil] 6583A ones, which are emitted in the out- 
skirts of the nebula, and sampling the largest variety 
of radial velocities. Note that the centered square 
2"x2" aperture misses the zone of zero radial veloc- 
ity, i.e. the zones where the velocity vector is entirely 
in the plane of the sky, while the rectangular aper- 
ture does not, leading to a very different line profile. 
In the case of the [O III] 5007A line, the rectangular 
slit actually intercepts a large zone with zero radial 
velocity, resulting in a profile without line-splitting. 
The H/3 line profile is the most featureless, essen- 
tially because of the strong effect of thermal velocity 
on this line. The important difference between the 
profiles of lines emitted by various ions amply justi- 
fies the use of a photoionization code to analyze the 
kinematics of real nebulae. 

The position of the slit with respect to the center 
of the nebula has a dramatic influence on the ob- 
served line profiles, as can be seen by comparing the 
dashed and solid curves in Fig. This means that 
a perfectly symmetric nebula can be mistaken for a 
nebula with uneven density distribution or complex 
velocity field simply because the slit is not perfectly 
centered! 

4. MIMICKING THE EFFECT OF A 
TURBULENT VELOCITY FIELD 

Re cently, iNeiner. Acker. Gesicki fc Szczer ba 
l)200fj) and iGesicki et all l(2003h have claimed to 
find evidence for turbulence in certain planetary 
nebulae, especially in those ionized by [WC] type 
central stars. While it is quite possible that the 
strong wind from [WC] stars induces turbulent 
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Fig. 2. Monochromatic images obtained for three emission lines (H/3, Hell 4689A, [Nil] 6583A, and [O III] 5007A), for 
the bipolar nebula shown in Fig. 0and described in Sec. 3.1. The apertures used to compute the emission line profiles 
shown in Fig.[3]are superimposed to the images: the sizes of the apertures are: l"x3" and 2"x2". The largest aperture 
(10"xl0") covers the entire nebula and is not shown. The centered apertures are drawn with solid lines. Off-center 
apertures are drawn with dashed lines. Surface brightness units are arbitrary, but the same for the four images. 
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Fig. 3. Some physical parameters, surface brightness distribution and line profiles for the bipolar nebula shown in Fig.Q 
and described in Sec. 3.1. Upper row: radial distribution of the hydrogen density (left panel), the electron temperature 
(middle panel) and the expansion velocity (right panel). Solid curves correspond to the polar direction, dotted curves 
to a perpendicular one. Second row : from left to right, surface brightness distribution in the H/3, [Nil] 6583 A and 
[O III] 5007A lines along the polar axis (solid curves) and in the perpendicular direction (dotted curves). Last three 
rows: from left to right, line profiles H/3, [N II] 6583A and [O III] 5007A lines, through apertures represented in Fig. |21 
(the size of the aperture is specified on top of each plot). Solid curves correspond to centered apertures, dashed curves 
to off-centered apertures. Intensity units are arbitrary, but the same for all the plots. 
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Fig. 4. Spherical nebula, with a Hubble flow expansion law and with negligible turbulence. The setup of the figure is 
the same as for rows 2 - 5 of Fig. [3] 
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Fig. 5. Spherical nebula, with a low velocity Hubble flow expansion law and high turbulence. The setup of the figure 
the same as for Fig. 2] 
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Fig. 6. Ellipsoidal nebula, seen pole-on, with a Hubble flow expansion law and without turbulence. The setup of the 
figure is the same as for Fig. 0] 
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motions in the nebular gas fsee lMellemall2003|) . we 
feel that the observational evidence is not as strong 
as might be thought. 

The effect of turbulence is to broaden the profiles 
by the same amount for lines of all elements. In ex- 
treme cases, this may result in the line splitting to 
disappear, if the broadening by turbulence is impor- 
tant. This is illustrated in Figs.0]and [51 which show 
the results from two spherical models. Both models 
are defined by an ionizing star of T o ff=70 kK and 
L*=l x 10 37 erg s -1 surrounded by a spherical nebula 
with an inner cavity of 7.5xl0 16 cm and a constant 
hydrogen density of 8000 cm -3 . The first model has 
a linear expansion law reaching a value of 40 km s _1 
at the recombination front, while the second one has 
a turbulence of 20 km s^ 1 added and the velocity at 
the recombination front is reduced to 25 km s _1 . In 
order to construct the surface brightness and veloc- 
ity profile figures, a distance to the object of 1.16 kpc 
has been assumed, implying an outer angular radius 
of 8 ". The setup of the figures is the same as for 
Fig-Olexcept that the first row has been omitted. We 
see in Fig.GDthat turbulence has strongly reduced the 
two peaks that were seen with the 3"x4" slit. 

There are however several other ways to achieve 
such a broadening. Even in spherical symmetry, one 
can produce similar profiles by changing at the same 
time the density distribution and the expansion ve- 
locity field. Another way is to consider an ellipsoidal 
nebula, rather than a spherical one, and oriented 
pole-on. With such a configuration, the weight of 
the zones of nearly zero radial velocity becomes large, 
and high velocity wings are produced as well. Fig.El 
shows such a case, which differs from the case shown 
in Fig. 21 in that the nebula is now ellipsoidal with a 
cavity of 7.5 x 10 16 cm and 1.12 x 10 17 cm and an inner 
density of 8000 cm" 3 and 3550 cm~ 3 , in the equato- 
rial and polar directions respectively. The expansion 
law is identical and reaches a maximum velocity at 
the polar ionization front of 70 km s _1 . As seen in 
Fig. [f)| this ellipsoidal model has very similar line 
profiles to the turbulent spherical model, at least if 
the slit size covers the whole nebula (and for H/3, 
even if the slit covers a substatial part of it and is 
off-center). Note that even the surface brightness 
profiles are very similar. Therefore, nebulae that 
would correspond to such models would be impos- 
sible to distinguish in practise, if using observational 
setups such as those represented by the last two rows 
of Figs. Eland El Note, however, that with a smaller 
slit, the two models can be distinguished, at least 
if the slit is perfectly centered (compare the second 
rows of Figs. 2]and HJ) because in that case the split- 



ting of the emission lines due to crossing two emitting 
regions on the line of sight. 

5. DISTINGUISHING A BLISTER FROM ITS 
SPHERICAL IMPOSTOR USING LINE 
PROFILES 

In Paper I we showed an example where two to- 
tally different geometries could lead to the same sur- 
face brightness H/3 maps: a face-on blister could be 
mistaken for a sphere - which we named the "spher- 
ical impostor" 3 . Here we illustrate how high resolu- 
tion spectroscopy can allow one to discriminate be- 
tween the two geometries. Indeed, the velocity fields 
for these two geometries are expected to be very dif- 
ferent. For the face-on blister, the gas is streaming 
out of the neutral cloud towards the observer, while 
for the spherical impostor, one expects a spherically 
symmetric expansion. 

Fig. shows the line profiles of the blister model 
of Paper I, in which we have assumed that the veloc- 
ity is perpendicular to the gas surface and is set to 
30 km s _1 . Fig.|Hlshows the line profiles of the spher- 
ical impostor model from Paper I, in which we have 
assumed a constant velocity of 30 km s _1 . The dif- 
ferences between the observed profiles between Fig.0 
and Fig. |H| are spectacular. For the blister, all the 
lines are blue-shifted and narrow, since gas flows in 
only one direction. For the spherical impostor, when 
using a thin slit passing through the center, the emis- 
sion lines are splitted. If the slit covers a large frac- 
tion of the object or if it is off center, there is no 
splitting, but the lines are broad (their width at the 
base is twice the expansion velocity). 

6. GENERATING 3D MAPS AND 
PV-DIAGRAMMS 

From examples shown in the previous sections, 
it is clear that the interpretation of line profiles is 
by no means straightforward. One of the problems 
is that observed profiles strongly depend on the po- 
sition and size of the slit. Obviously, a large amount 
of information is necessary to unambiguously unveil 
the morphology and kinematics of a nebula. There 
are now an increasing number of instruments which 
use integral field units (IFU) to achieve spectroscopy 
of an extended portion of the sky, either using lens 
arrays, optical fibres or image slicers. These in- 
struments can thus generate 3D maps of extended 
objects, with wavelength (or velocity) being the 

3 In the right panel of Fig. 17 in Paper I, the line styles 
to represent the blister and the spherical models have been 
inverted: it is the blister model which was represented by a 
continuous line, contrary to what is stated in the legend and 
to what has been done for the other panels of the figure 
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Fig. 7. Line profiles for the face-on blister (see Sect.|K| The setup of the figure is the same asfor rows 3 - 5 of Fig. |S| 




Fig. 8. Line profiles for the spherical impostor (see Sect.^J. The setup of the figure is the same as for rows 3 - 5 of 
Fig. El 
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Fig. 9. A 3D map of the bipolar nebula model from Sec. 13.11 Each quadrant corresponds to a different emission line. 
The surface brightness is represented by levels of grey. Superimposed on this image are the profiles of the same lines, 
integrated over areas of 5x5 pixels. 
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third dimension. Such examples can be found in line by assuming an arb itrary velocity field. In this 

lAmbrocio-Cruz. Laval. Rosado. Georgelin. Marcelin. Com OTOik.Bhlrrotte Asc¥iAte g enerate the photoio nization 
(20041: IVasconcelos. Ceraueira. Plana. R.aga fc Morissed model was NEBU_3D IjMorisset et al.ll2005(l but any 
l|2005h . Our tool is well-suited to also produce such 
maps for photoionization models of asymmetric 
nebulae. One such example is shown in Fig. 
which represents the bipolar nebula model from 
Sec. 13.11 Each quadrant corresponds to a different 
emission line (H/3, Hell 4689A, [Nil] 6583A, and 
[OIII] 5007 A). The grey-colour image represents 
the surface brightness in the line, computed with 
pixels of 0.075"x 0.075". Superimposed on this 
image are the profiles of the same lines, integrated 
over areas of 5x5 pixels. The line intensities arc all 
on the same scale. Apart from illustrating how the 
line profiles depend on the emitting ion, this figure 
also clearly shows how the profiles change with 
position. In the present case, they may change from 
single-peak to triple-peak structures. The latter 
occur near the waist of the bipolar structure, when, 
due to inclination, the line of sight crosses the front 
and back of one bubble, plus the front on the other 
bubble (see Fig. \I\for visualization). Note that the 
largest component of the double-peak line profiles is 
the red one in the upper spectra and the blue one in 
the lower spectra because of the inclination of the 
nebula. 

To illustrate the capacities of VELNEB_3D to 
produce Position- Velocity (PV) diagrams, the bipo- 
lar model described in Sec. 13.11 is used to generate 
images presented in Fig. ^] The H/3 surface bright- 
ness map is shown with the two narrow slits superim- 
posed, which are used to compute the PV-diagrams. 
One is obtained for an horizontal slit, the second 
one for a vertical one. They are both centered, but 
any position can be defined by the user. The nebula 
has its polar axis making an angle of 45degr with the 
plane of the sky, resulting in well defined red and blue 
components for the vertical slit (upper right panel of 
Fig. I10|) . when the slit is crossing the two bubbles 
of the nebula. This orientation effect is not seen in 
the case of the horizontal slit crossing mainly gas at 
projected velocity closed to zero (lower left panel of 
Fig. Iiop . The last panel shows a velocity channel 
map obtained for velocity close to 11 km s -1 . The 
part of the nebula expanding in the positive direction 
is well traced. 

7. SUMMARY AND PROSPECTS 

In this paper, we have presented a tool, VEL- 
NEB_3D, which uses the results from any 3D pho- 
toionization model to generate emission line profiles, 
position-velocity maps and 3D maps in any emission 



other code can be used. 

So far, most of of the modelling of line profiles 
in p lanetary nebulae h as been done using ID codes 
(see lGesicki et al.l2003l and references therein). The 
examples given in this paper show that by taking 
into account deviation from spherical symmetry, the 
interpretation on the velocity field may be different. 
This is important, s ince about 80 % of planetary neb- 
ulae are not round (Balick 2006). 

With this tool, we have been able to show how 
much the interpretation of observed line profiles may 
depend on the exact position of the slit. For exam- 
ple, a complex line profile may be obtained even with 
a simple expansion law if the nebula is not spheri- 
cal and the slit slightly off-center. Since physically, 
geometry and velocity field are related (the geome- 
try being the result of the evolution of the velocity 
field in time) it is likely that in real nebulae with 
complex profiles, both the geometry and the veloc- 
ity field are in fact complex. Our tool allows one 
to explore a larger parameter space than ID models 
and (perhaps) pin down the best solution for ob- 
served nebulae. It can also be useful as a tutorial 
kit for those wishing to better understand observed 
monochromatic images and line profiles. 

We believe that trying to reproduce obser- 
vations of real ionized nebulae with our tool 
should give more insight into the physics of these 
objects and guide the evolutionary dynamical 
modelling of planetary nebulae such as done 
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While dynamical modelling is obviously the ulti- 
mate step in understanding the origin and evolution 
of planetary nebula shapes, or of the shapes of other 
kinds of nebulae, it involves heavy computation and 
at the same time a simplification of the physics at 
work. Therefore an exploration with a tool such as 
we propose is worthwhile. 

Our examples suggest that it may be actually 
very difficult to reconstruct the geometry and veloc- 
ity field of any real nebula without any underlying 
assumptions. However, there are cases where high 
resolution spectroscopy may distinguish between two 
very different geometries. Such is the case of a face- 
on blister and its spherical impostor, which have the 
same surface brightness distribution in Ha for the 
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Fig. 10. PV diagrams obtained for the bipolar nebula presented in Sec. l3.1I IJpper left: surface brightness image for H/3 
with the 2 slits used for the forthcoming PV diagrams, Upper right: PV diagram obtained for a narrow vertical slit 
crossing the center of the nebula, Lower left: PV diagram obtained for a narrow horizontal slit crossing the center of 
the nebula, Lower right: velocity channel map obtained from the V=ll km s _1 channel. 
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same ionizing star. 

Our tool may also be useful in preparing and in- 
terpreting observations of ionized nebulae with inte- 
gral field units, which are able to provide 3D maps 
and will predictably lead to a major progress in neb- 
ular astronomy in the coming years. We foresee 
however the difficulty that with such a huge amount 
of observational and computational data, it may be 
tough to pin down the really important physics, in 
other words, to not only answer the question "what 
is it that we see?" but also "why is it like that?" 
Clearly, ways will have to be found to make the best 
use of these tools. 

Another application of our tool is the possibility 
to compute a large amount of models, varying the 
morphology, the ionization parameter, the velocity 
law and the position of the slit and to include the 
resulting emission line profiles in a Virtual Observa- 
tory (VO). The users of the VO would be able to 
quickly scan the catalog to search for all the nebulae 
satisfying given criteria, for example on line profiles. 
Such a work is in progress. 

The VELNEB_3D tool is available on request to 
CM. 
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UNAM, Mexico, for hospitality and financial sup- 
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